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KN interactions are investigated via an effective non-linear chiral meson-baryon Lagrangian. The 
adjustable parameters are determined by a fitting procedure on the K~p threshold branching ra- 
tios and total cross-section data for p l £ b < 250 MeV/c. We produce predictions for the mass 
spectrum, and scattering lenghts a K - p , a n (K~n — > K~n), a„(K n —>K n), and a ex (K~ p —tK n). 
The KN amplitudes thus obtained, as well as those for other two-body channels (ttN, NN, and 
YN) are used as input to predict the scattering length A K - d , for which we have devised a rela- 
tivistic version of the three-body Faddeev equations. Results for all two- and three-body coupled 
channels are reported both in isospin and particle bases. All available KN data are well reproduced 
, and our best results for the K~p and K~d scattering lenghts are a K - p = (—0.90 + i0.87) fm and 

(N ■ A K - d = (-1.80 + il.55) fm. 

PACS numbers: PACS numbers: 11.80.-m, ll.80.Jy, 13.75.-n, 13.75.Jz 
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I. INTRODUCTION 



The present paper is devoted to the study of the -fT~-deuteron scattering length A^-^ by exploiting a relativistic 

j» version of the three-body Faddeev equations in which the principal two-body input is based on an effective non-linear 

chiral meson-baryon Lagrangian in the strangeness S = — 1 sector. A preliminary version has been reported in Ref.pJ. 

■ Below we shall begin with a survey on the two-body KN amplitudes which are the central input to our present 
enterprise, as well as the evolution of the low energy K~ d physics. In this way we hope that our motivation be clear 

■ to the reader. 

While the low to medium energy kaon-nucleon processes (say, pj^ ab < 1.5 GeV/c) have been known to show no 
significant structure up to the pion production threshold in the KN channel, its u-channel counterpart: the KN 
channel, presents quite a rich structure (resonances, possible bound states in the continuum, etc.): see for example 
For the most recent discussions on both experimental and theoretical fronts, see, for example, Olin and Park Q. 
' The thus mentioned characteristics may be understood within a simple quark model where the quark structure of K, 
K and N are known to be qs, qs, and qqq, respectively with q being u and/or d quarks. In this picture iheKN system 
may be rearranged to become a combination like {qq)(sqq) which may be identified, for example, as 7rY, (Y — A, E) 
in terms of the lowest-lying octet hadrons, or to generate strangeness S = — 1 hyperon resonances. On the contrary, 
such a scenario does not materialize for the KN system in which the antiquark involved is s, hence the corresponding 
low energy process is uniquely the elastic scattering of K and N. So both theoretically and experimentally, the KN 
(S = —1) system has been drawing much more attention than the KN (S — 1) channel. 

Of particular interest in this regard has been the K~p channel near threshold. It is dominated by the below- 
threshold resonance A(1405) to which it strongly couples. This resonance decays almost exclusively to 7rE. One of 
the intriguing subjects related to this resonance has been on its dynamical origin: whether it is a hadronic bound 
state of K~p embedded in continuum (since it is located above the irY, threshold), or a three-quark baryon resonance 
(or something more exotic) . Whereas no definite conclusion has been drawn based on the scattering data analyses by 
forward dispersion relation 2, 3] , a recent effective chiral Lagrangian approach Q has presented a convincing picture 
in favour of the K~p (I = 0) bound state, as we shall discuss later in the context of the objective of the present paper. 
A more down to the earth, but quite important problem has been the K~p scattering length o^- p , that dictates the 
threshold characteristics of the K~p interaction. For this quantity the so-called kaonic hydrogen puzzle, see e.g. Ref. 
0, had disturbed the community engaged in low energy meson-baryon interactions for quite a long time. Briefly, 
the puzzle originated from the fact that the real part of a^-- p extracted from the Is atomic level shift (due to the 
strong interaction) of the kaonic hydrogen had an opposite sign to the one obtained from the analyses of the K~p 
scattering amplitude. 



Despite the persistence of this puzzle, several pioneering works on the low energy negative kaon-deuteron (K~d) 
scattering by solving three-body equations were performed by simply disregarding the information from the kaonic 
hydrogen. The first study of that type on the K~d elastic scattering at low energies was performed as early as 1965 
[8j by adopting simple iS-wave rank one (non-local) separable interactions for the 7 = 0,1 KN channels and for the 
3 Si deuteron channel as basic two-body ingredients. For the KN interaction the potentials were fitted to reproduce 
the complex- valued K~p scattering length known from the amplitude analysis at that time (see for example Q). 
The principal objective of the work was to see the convergence properties of the multiple scattering series: at low 
energies they found that the single and double scattering contributions were far from sufficient. An extension of this 

work to the three-body break up channels: K~d — > K~np, A nn was performed by the same authors by deforming 
the momentum i nteg ration in the complex plane to avoid the anomalous threshold branch cuts due to the final state 
A A interactions [l(J- The result turned out to be insufficient for discriminating the different K~p input amplitudes 
by comparing with the data, as the two-body input were incomplete, along with insufficient statistics for the data. 
However, the methods developed served in motivating later attempts for studying the three-body final states from the 
theoretical side. Also the same authors improved their first elastic scattering calculation in Ref. |8( by incorporating 

the isospin breaking effect manifested in the mass difference between A° and K~ as well as between n and p, which 
resulted up to an about 10% difference in the cross section Later Myhrer ^3] studied the role of the A(1405) 
resonance in the threshold K~d scattering by assuming a simple resonance form for the input K~p amplitude. Here 
again the insufficiency of few first iterations of the multiple scattering series was found in the presence of a two-body 
resonance. This was then followed by using improved A N amplitudes that took into account the effect of coupling to 
the physically accessible irY, (Y — A, E) channels and, though implicitly, the effect of A(1405) 13]. In this regard one 
should be reminded that in the earliest works mentioned above ja, EI those effects were implicitly represented only 
by the complex-valued KN scattering lengths. To summarize, all the models mentioned here summed up multiple 
scattering series driven by an S-wave K N scattering with a spectator N and a 3 <5i NN scattering in the presence of 
a AT spectator. We will refer to this type of models as Sing le- Channel Approach. 

The next generation of theoretical endeavour [lj, [jjj, \wL Il7| may be characterized by explicitly taking into account 
the three-body channels involving hyperons: irNY, (Y = A, E), where Y is produced from the A A — > irY reactions 
in the presence of a spectator nucleon. Also included were the two-body interactions nN with a spectator Y as well 
as YN with ir as a spectator. In this way genuine three-body unitarity was guaranteed to hold. Common to all three 
works cited here is the way to construct all the two-body input amplitudes to the three-body equations. Apart from 
the NN deuteron channel, all the two-body amplitudes: coupled A N- irY channels, irN, as well as coupled NY-NY 
channels, were assumed to be obtained from rank one separable potentials (mostly in S'-waves) where the strengths 
and ranges were adjusted to fit the available cross sections, etc. Isospin was assumed to be conserved exactly so that 
the number of channels to deal with in the two-body input as well as in the three-body equations be contained as 
manageable. This will be termed as the Multi-Channel Approach as compared with the single-channel one mentioned 
earlier. 

The first two of those multi-channel approaches [T3. ITHj l were dedicated to the near threshold break-up reactions: 
K~d — > nNY, where data exist for the reaction rates, the final mass spectra m(pA), the neutron momentum spectra in 
the final nE~7r + as well as in nE + 7r~ three-body channels [T^ . [l9| . The models were able to reproduce the experimental 
trends depending upon which combinations of signs in the non-diagonal amplitudes such as A A — > 7rS, A A — > EA 
should be adopted. Note that those amplitudes determined by fit to the corresponding cross sections are unique up 
to an overall sign, or phase in more general terms, unless additional constraints from, say, some symmetries, were 
imposed; but in three-body processes the sign difference does show up! Also difference in the two-body amplitudes 
responsible for the final state interaction was found to be visible in the break-up channels, so the data could be used 
to discriminate different two-body models used. In both works the authors also calculated the K~d scattering length 
A K ~ d . The best value obtained for this quantity in those works may be identified as (-1.34 + i 1.04) fm 

The third multi-channel three-body calculation was carried out in Refs. [T^ . lr^ | for low energy elastic K~d scattering 
as well as to find the best theoretical value for A K - d . The data used to the x 2_n ts to determine the two-body separable 
interactions were the same as in | 15| . The novelty was to use relativistic formalism, hence the correct kinematics was 
ensured when dealing with different total masses (or different thresholds) in the entrance and exit channels, both 
in two-body input amplitudes as well as in three-body coupled equations. With the 3 Z?i partial wave component 
included in the deuteron channel interaction, the best value was A K - d = (—1.51 + i 1.45) fm. 

After the last theoretical calculation just discussed above had come out in 1990 ^(|, the low energy Kd physics 
became dormant for about ten years. One may identify some of the possible reasons for that void: 

(i) The kaonic hydrogen puzzle mentioned earlier kept persisting, so it was felt that without any solution to it, one 
could not find any credible low energy K~p amplitude for use to improve K~d models. 

(ii) The rank one separable potentials adopted to model the essential ingredients, namely the coupled KN-iiY 
amplitudes, lacked support from underlying strong interaction theory. So, even by adopting isospin symmetry to fit 
the existing data, there was no compelling reason to believe that the best fitting amplitudes be really acceptable on 



physical ground (a serious attempt of this type dated back to Henley et al. and Fink et al. 1201 '). In this respect some 
efforts to constrain the fitting parameters by SU(3) symmetry deserve to be noted pll l22j. However, the separable 
ansatz still needed to be given a proper justification. It should be useful to note that there were also local versions 
of the corresponding potentials based upon physical constants from Chiral Lagrangian to be somewhat tuned [23| . 
but the local form had no support from the supposedly more fundamental theory, either. In addition to this, taking 
into account the isospin breaking effects would have pushed the picture farther into the mist. So why should one go 
forward under such circumstances? 

Quite fortunately, there were two major breakthroughs: one on the experimental and the other on the theoretical 
front. 

First, the long haunting kaonic hydrogen puzzle was finally put to an end by the KEK experiment |24| . and the K~p 
scattering length a c K _ p , where Coulomb effects are not separated, was extracted to be: a x=jL = ( — 0.78± 0.l5(stat) ± 
0.03(syst)) + i(0.49± 0.25(stat) ± 0.12(syst)) fm, by applying the Deser-Trueman formula |25j to the kaonic hydrogen 
Is level shift T, and width e, by K~p strong interaction : 



e + i-=2a 3 fi 2 a c K - p , 

with a the fine structure constant and [i the K~p reduced mass. Although the thus obtained quantity includes 
the effect of the Coulomb interaction, hence not identical to the corresponding quantity due exclusively to strong 
interaction, the difference is, even conservatively, at most within a few percent: in the case of pionic hydrogen the 
difference appears to be below 1%, see for example [2(| I21J, so the extracted He{a, K- v ) finally was found to have 
the same sign as the one from the analysis of scattering data. The DEAR project [23] with the DAQNE facility 
(0-factory) at Frascati had been planning to repeat the experiment with higher precisions, and the data taking is 
reported to have been over. In addition, along with other experiments involving various if-mesons, this project has 
planned to measure the corresponding quantities for the K~d atom, and the experiment is expected to start soon 
|29j . To obtain the K~n scattering length without recourse to isospin symmetry (hence to find out how good that 
symmetry is realized in the kaon sector in this quantity) is one of the objectives for this measurement. But, a more 
ambitious picture such as using the quantities obtained to extract the kaon-nucleon Sigma-term: o~kn, with the help 
of theories such as Chiral Perturbation Theory, see for example Ref. [30j, has been some strong driving force for the 
experiment E l l8ll l32j. In fact, such a program has been put into practice recently using the pionic hydrogen and 
deuterium [23j. It should be useful to mention here that extracting the Coulomb corrected K~d scattering length 
through the Dcscr-Trucman formula may not be as accurate as that for the K~p system, and that relation between 
the purely strong and Coulomb corrected K~d scattering length may not be quite simple either. This was discussed 
in terms of simple models by Barrett and Deloff |2J|, see also [2^. This will be revisited in Section Ivl . 

Second, based upon chiral perturbation theory, there has been a steady progress in describing the low to medium 
energy meson-baryon interaction. And within the context of our present interest, there was an important breakthrough 
made by Oset and Ramos in describing the coupled meson-baryon channel amplitudes for the S = — 1 sector |5(- The 
channels involved are: KN, irY, and rjY. The driving (potential) terms to the two-body coupled Bethe-Salpeter 
equations were taken from the lowest order in the effective chiral Lagrangian for the 0~ octet mesons and l/2 + octet 
baryons. Then, an on-shell ansatz developed in |35j was introduced, which enabled the authors to transform the 
coupled integral equations into a set of algebraic equations where the major task was to deal with the integration of 
various meson-baryon two-particle propagators which are ultra-violet divergent. With only two parameters adjusted 
to very plausible values, the cut-off at p = 630 MeV/c in the momentum integrations to be convergent, and the 
average octet meson decay constant set as / = 1.15/^, the authors were able to reproduce various low energy data 
associated with the above-mentioned coupled channels impressively (the isospin breaking effect in the particle masses 
was included). Particularly, the A(1405) was generated dynamically as an unstable bound state of K~p in the 1 = 
channel whose decay into ttY was also correctly reproduced. As will be discussed in the next Section, what made 
this approach particularly noteworthy is that the on-shell ansatz they adopted turned out to be interpreted as a 
practical justification for a separable representation of two-body meson-baryon potentials, where coupling strengths 
are the product of two parts, one dictated by SU(3) symmetry, the other depending on energy. Note that this on-shell 
ansatz was made to be rewritten in a more elegant form, viz. in the so-called N/D representation, or cast into a 
once-subtracted dispersion integral representation of the two-particle propagators, etc. There, the dynamical left-hand 
cut contributions were shown to be weak, and the divergent integrals were made finite by an introduction of several 
subtraction constants adjusted to reproduce the relevant data, including the KE channel [Sfil Isfl IssL 13^. liflT] . 

Now, upon witnessing the experimental and theoretical progress reviewed above, the time is ripe for starting the 
K~d scattering study again. In fact, the Oset-Ramos result y| was applied to calculate the scattering length A K - d 
within what is called the fixed center (or fixed scatterer) approximation (FCA) by Kamalov et al. |4lj . Here we want 
to claim that a more refined approach should match the Oset-Ramos type amplitudes. More specifically, we think it 



necessary to exploit a reliable method in dealing with three-body scattering. 

In this article, we thus present a complete study of the K~d scattering length within the Faddeev equations. In 
Sec. Ill Al we first adopt the approach due to reference [j| and construct the coupled s-wave separable -ftT TV interactions, 
including the coupling to the charge exchange, as well as the ttY and r/Y channels. We study different models, with 
parameters fitted to the available experimental data, along with the constraint to remain close to the SU(3) values. 
The results are presented both in the isospin and particle bases where the latter takes into account the isospin breaking 
effect in terms of the physical meson and baryon masses. The parametrizations chosen for the deuteron channel are 
described in Sec. Ill Bl Then we discuss other two-body interactions adopted in our study in Sec. Ill CI In Section ITTT1 
we review the structure of the relativistic three-body Faddeev equations. In Sec. IIVI we calculate the K~d scattering 
length by solving those equations, both in isospin and particle bases. To our knowledge, this is the first calculation 
done in the particle basis, permitting to evaluate the isospin breaking effects at the three-body level. The sensitivity 
of A K - d to the two-body input is investigated. The discussion is developed in Sec. [V] and our conclusions are given 
in the last Section. 



II. TWO-BODY INTERACTIONS 

In this part, we describe the various two-body interactions used as input in the three-body equations. In the 
first two Subsections, we present the K N and NN interactions, which are the fundamental ingredients to the K~d 
problem. The last Subsection is devoted to a brief description of the irN and YN interactions. 

A. EN INTERACTIONS 

Here, we intend to construct two-body coupled channel AT AT interactions. The coupled channels involved are K N, 
ttY, and rfY (Y — A, S) with different total charge states for the mesons and baryons (in the particle basis), or in 
different total isospin states (in the isospin basis). See Appendix lAlfor the relation between those two representations. 
The physical masses used in the particle basis and average masses used in the isospin basis may be found in Table [i] 



TABLE I: Particle masses (in MeV). The fourth line gives the average mass for each isospin multiplet, and the last line specifies 
the phase convention used for the isospin states. 
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1. Separable models 

Let us use i, j, k..., etc. as channel indices. Since our present interest is in the K p near its threshold, we may 
safely assume that any given meson-baryon system in the coupled channel is in the relative orbital angular momentum 
s-state [l = 0). So we may adopt the s-wave projected coupled Bethe-Salpeter (or relativistic Lippmann-Schwinger) 
i-matrix equations |4^ | for the transition j — > i which takes the following form: 

T ij = Vi J +J2VikG^T kj , (1) 
k 

where Vij is the transition potential, and Gg is the free meson-baryon propagator for the intermediate channel fc. We 
note here that implicit in the above expression are that (i) the meson-baryon systems are in the center of mass frame, 
and (ii) the integration is performed over the off-shell four momentum associated with channel fc. 

We take two additional simplifications to make the coupled equations manageable. The first one is to adopt the 
Blankenbecler-Sugar procedure to reduce the momentum integration from four to three dimensions |43l l44l |45| . In 



particular, the two-particle propagator is re-expressed as G§ = GqIj^ct), where pk is the magnitude of the three- 
dimensional relative momentum of the intermediate channel k, and a is the square of the total center of mass energy. 
This may be done by taking the discontinuity of Gq over the unitarity branch cut and use it to represent G§ in a 
dispersion integral form. We note that this procedure results in on mass shell but off energy shell form of equations 
which may be regarded as a relativistic extension from a familiar non-relativistic scattering theory. The second 
simplification step is to assume that the s-wave potentials take a non-local separable form: 



^^SiW^Wfljfe), / = 0,1 , (2) 

where I is the total isospin for the meson-baryon system, and gt is the cut-off form factor for channel i which is 
assumed to be a function of the magnitude of the three dimensional relative momentum vector in the same channel. 
In general the coupling strength Xjj is assumed to be a function of a as indicated in the above equation with no 
left-hand cut assumed. Some rudiments of how our coupled two-body i-matrices may be obtained with the separable 
interactions is found in Appendix iBl 

Now, when we compare the expression for the coupled KN channel t-matrices, Eq. (20) in 5] with our coupled 
channel i-matrices, Eq. I)B6|1 for which the related quantities are defined in Eqs. (|B7I IB8JI . we see immediately 
that the two results are identical provided that (i) we set gi(pi) = 1 for all i, (ii) impose a momentum cut-off p ma x in 
the integration in Eq. (|B8fl . and (iii) set 



^•=-^475^ + ^)- (3) 

The expression for X\a above is from Oset and Ramos [5j, which was obtained from the lowest order expansion in 1/ / 
of the chiral Lagrangian for the octet 0~ mesons coupled to the octet l/2 + baryons. The coefficients C\a are due to 
SU(3) symmetry and tabulated as: C^ =0 = £)y, Cjp 1 — Fij in that publication. These are convenient for the isospin 

basis, but may be trivially transformed to the corresponding coefficients CV/ (for K~p) and Cy (for the K~n related 
channels), for use in particle basis, also tabulated in [5j. The corresponding change to obtain the strength parameters 
in the particle basis in terms of the Xj- follows trivially, see for example Eq. (|A2Jl . In our equation ® above, ej and 
€j are the meson energies in the centre of mass system for the i and j channels, respectively. 

Though our argument above went just in the opposite direction to what one finds in Oset and Ramos [f| (see also a 
more formally trimmed version of the Oset-Ramos line of derivation by Nieves and Arriola H^), we have established a 
practical equivalence between the separable potential and the on-shell ansatz for the coupled meson-baryon equations. 
So within the framework of effective meson-bayon field theory, it is now possible to claim that the separable ansatz is a 
very reasonable starting point in describing the s-wave interactions at low energies. In our present work we choose to 
retain the form factors rather than imposing a sharp cut-off. This is due to the fact that when solving the three-body 
equations we rotate the momentum integration path off the real axis and into the complex plane. For that purpose, 
a sharp cut off is not practical. Then in order to respect SU(3>) symmetry, we choose to use a single form factor for 
all the different channels with a monopole form: 



2 

9(P) = Z2X*2* (4) 
p z + p 

where (3 is the effective cut-off momentum. Based upon the discussion above we adopt two slightly different types 
of interactions. The first one which we call OS1 is just V^, Eq. (J2J), with = gj = g(p), and Xjj from Eq. 
We expect this interaction to produce a very similar result to the original Oset and Ramos model. The second 
interaction model is called OSA, which is a variant of OS1 in that it incorporates the possible SU(3) breaking effect 
in the coupling strengths (or in the meson decay constant /) in terms of extra parameters with the substitution: 
C\a — > b\jC\y We then adjust on the relevant data f3, /, and &L's, the last ones only for OSA, with a constraint that 
the SU (3) breaking effect is reasonably contained. 

It is important to stress here that we follow the observation by Oset and Ramos Q and retain the r]Y channels in 
our fit although this channel has a substantially higher threshold as compared with that for KN, see Table [I] The 
necessity for the inclusion of these channels will be demonstrated later. As a result we have three coupled channels: 
KN, 7rE, 77A to deal with for I — 0, and four coupled channels: KN, 7rE, it A, 77S for / = 1. In terms of physical (or 
particle) channels the following two groups are separately coupled: 



K~p , K n , Att° , S H 



SO J 
7T 



, At? , I] /?, 



(5) 



R-n -> R-n , Att" , £°7T- , ITvr , E~r?. (6) 

Below in the next Subsection we briefly review available experimental data to which the model interactions are fitted 
and compared. 

2. Experimental data 

Since there are no data associated with the initial K~n channels in Eq. ©, we only discuss the ones in Eq. JSJl. 
Furthermore, in the low energy range of our current interest, i.e. p l £ b < 250 MeV, the last two channels involving 
the r\ meson are physically closed as their thresholds are substantially higher than the rest. Otherwise the remaining 
physically accessible coupled channels are now strongly influenced by the / — A(1405) resonance below the K~p 
threshold which decays almost exclusively to 7rE. Note also that while the K°n has a slightly higher threshold, all the 
irY channels have lower threshold than that for K~p. So altogether there is a very rich structure in this coupled set of 
channels 1 .. In the energy range considered here, some 90 data points are available EE E3, El, E^, H3, E3, [HE E3, E3 
These data, obtained between 1965 and 1983, bear unequal accuracies, as briefly discussed later. 

Moreover, accurate data EE EE EE EE E3 f° r threshold branching ratios are also availabe, i.e., 

7= *n 3^^ = 2-36 ±0.04, (7) 

fc^O 1 (A p — > 7T E+) 

T(K~p — > charged particles) . . 

R c = lim l /T / -| f r-i =0.664 ±0.011, 8 

fe^o T(K-p — » all final states) 

Rn = lim— F(X I — n ° A } r = 0.189 ±0.015. (9) 

fe^o T\K p — > all neutral states) 

There are also data 58] on the invariant mass spectrum of the E + 7r _ system, which have been exploited to investigate 
the nature of the A(1405) resonance. 

Finally, the last piece of crucial experimental information comes from the recent KEK measurement [24| of the 
K~p scattering length. As explained in the Introduction, the obtained value which includes the Coulomb effect: 
a c K - = (-0.78 ± 0.15 ± 0.03) + i(0.49 ± 0.25 ± 0.12) fm, resolves the "kaonic hydrogen puzzle". 



3. Results of the fit 

For model OS1 interaction, we have adopted the same strategy as in Ref. jE| by fitting our parameters / and j3 
to the threshold branching ratios, with / constrained to deviate from / ff by less than ±20%. All other observables 
are "predicted", i.e. they are evaluated with the values of the parameters reached at the end of minimization (these 
values can be found in Ref. [lj). The branching ratios and the ax- p scattering length obtained in this model are 
comparable to the values from the Oset-Ramos model, see Table ITT1 

The same conclusion holds when we compare the total cross sections and the 7rE mass spectrum given in Figs. ^ 
and [3 with the correponding results in Ref. ||- 

For the OSA interaction, the parameters are /, (3, and bjj's, which are fitted to the threshold branching ratios and 
the K~p initiated cross sections, with the bfj's constrained to remain within ±30% of their exact SU(3)-symmetry 
values. The obtained values are given in Table ITTT1 2 .. 

Note that in the lowest order chiral Lagrangian approach, some coupling coefficients: C^ = D^j, Cjp 1 = as 
found in Tables 2 and 3 in Ref. [5j are equal to zero. We have chosen to keep these zero values, thus we do not need 
the corresponding SU (3) breaking coefficients 6f , which is reflected in Table II I II 



1 Notice that the contributions from the K~p K°S°, K + S~ channels, which couple only indirectly to the KN system, were found 
insignificant 0] at low energies, and hence are not included in our formalism. Moreover, our aim in this work being to use the ele men tar y 
reactions as input in the K~d elastic scattering, the electromagnetic processes K~ p — » A7 , S°7, studied by other authors |21, 22], 
are not considered here. 

2 In our previous publication [lj, interaction OS2 was introduced which took an identical form to the current OSA. The resulting fits by 
these two models are very similar. The major difference is that the 5(7(3) breaking has come out less in OSA than in OS2, due to a 
more careful handling of the data-base. 



TABLE II: K p threshold strong branching ratios and K p scattering length (in fm), calculated in the particle basis. The 
results in the isospin basis are shown in italic letters. See text for experimental data References. 



Authors [Ref.] 


7 


Rc 


Rn 




Im(a K - p ) 


Model 


Present work 


2.35 


0.651 


0.189 


-0.90 


0.87 


OSA 




3.17 


0.650 


0.257 


- 0.75 


1.11 


isospin basis 




1.04 


0.655 


0.130 


-0.74 


1.41 


OSA, r]Y excluded 




2.36 


0.657 


0.193 


-0.98 


0.80 


OSB 


Bahaoui et al.^] 


Z.oo 


U.DoD 


n 1 71 


— l.U4t 


U.00 


Uol 




O.OI 


U. VfCV 


n 9 / J 

U. <o>£4 


O'i 

— U . ao 


1 OR 


LSUSpill ULLSLS 


Oset & Ramos \£j 


2.32 


0.627 


0.213 


-1.00 


0.94 


Chiral rjY included 




3.29 


0.617 


0.292 


-0.85 


1.84 


isospin basis 




1.04 


0.637 


0.158 
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P, 1 " (MeV/c) P, 1 "" (MeV/c) 

FIG. 1: Total cross sections initiated by K~p, calculated in the particle basis with different K N models: OSA (bold full line), 
OS1 (bold dashed line), OSA-?; excluded (reg ular dotted line). The reg ular full line is obtained with model OSA in the isospin 
basis. Experimental data are from Refs. prtlH lisl El 15(1 Isl. I5I ImLIsII . 

Our results for OSA are presented in Table [H] and Figs. ^ 121 along with available data as well as values from a 
few earlier models. We first discuss the result for K~p initiated channels in the particle basis. As shown in Table ITT1 
the threshold reaction ratios: 7, R c , R n are better reproduced by OSA than by OS1. Regarding the K~p scattering 
length, our result is a prediction since we have not used the value extracted from the kaonic hydrogen atom |24| as 
part of the constraint in the fitting procedure. The real part is obtained closer to the data [24| by OSA than by 
OS1. As for the imaginary part, both models give the values at the limit of the experimental uncertainties within 2 
standard deviation. We point out that this trend for Im(aK- P ) is systematically observed with either the present 
separable models or with the Oset-Ramos approach. The cross sections are also better reproduced with OSA than 




FIG. 2: 7rE mass spectrum, obtained with the same models as in Fig. □ The bold dotted line is model OSB. The OSA results 
in the isospin and particle bases are practically identical. Experimental data are from Ref. 58]. 



TABLE III: Minimization results for the adjustable SU(3)-symmetry breaking coefficients of models OSA and OSB for the KN 
interactions. For OS1, these coefficients are equal to 1, the average meson decay constant is / = 1.20/^-, and the range is 870 
MeV. We use /,r=93 MeV/c. The reduced \ 2 are roughly 1.2. 



Isospin Channel OSA OSB 

/ 1.20 U 1.12 U 

KN 0.994 1.048 

KN-tyT, 1.108 0.917 

KN-rjA 0.851 0. 

ttE 0.903 0.926 

1 KN 1.056 0.833 

KN-tvZ 1.293 1.209 

KN-nA 0.943 0.933 

ttE 0.991 1.290 

KN-rjZ 0.757 0. 

range (MeV) 888.7 879.6 



with OS1, see Fig. ^ especially those for the K~p — ► 7r~£ + and the K~p —> K n channels which now have the 
correct magnitude as compared with the data. The position of the A(1405) resonance predicted by the present models 
are quite similar, and in good agreement with the data, see Fig. El (and also Fig. |1J. Finally, we have calculated the 
cross section of the K~p — > r/A reaction near threshold. We have found that the results predicted by models OSA and 
OS1 are in reasonable agreement with the recent data |5^] from the Crystal Ball detector at Brookhaven: the steep 
rise juste above the threshold, and the s-wave behavior expected for pseudoscalar meson production at threshold, are 
correctly reproduced. 



As for the results obtained in the isospin basis we only discuss the case with the OSA interaction, as the characteristic 
feature is the same with the result from the OS1 interaction. Here, we use the average masses of the hadron isospin 
multiplets everywhere (see Table QJ, hence now, for example, the charge exchange reaction: K~p — > K n becomes 
elastic energy-wise. Note that to calculate reaction cross sections, the entrance channel energies are taken using the 
physical particle masses, which are also used to find the final state phase space volumes. In Table ITT1 we observe large 
differences between the values of the threshold ratios 7 and R n calculated in the two bases. The shift of the cusps in 
the cross sections when using the isospin basis are clearly exhibited in Fig. On the contrary, the position of the 
A(1405) resonance is only slightly affected. These differences between the results in the two bases are similar to those 
observed in Ref. 0. 

Next we study the effect due to the contribution from the rjY channels. Given that the difference between the two 
thresholds for r/A and K~p final states is as large as ~ 230 MeV (with 77 £ the difference is about 75 MeV larger), 
hence a priori the effect from the coupling to these channels is expected to be insignificant. However, in Ref. |5j the 
authors observed the importance of retaining these channels. So we want to check their claim. Very qualitatively, 
the role of these channels may be best understood in the exact isospin symmetry limit. Then rjY channels are the 
only ones whose thresholds are above the one for the KN. Hence, they provide a definite attraction to the elastic 
KN process. As a result the coupling to rjY states controls the binding properties of the KN (in the effective chiral 
interaction adopted here, the 1 = A(1405) is a bound state oiKN embedded in the continuum state of the 7r£ 
channel). As a first step, we have discarded these channels from model OSA by forcing the KN-rjA and KN-rjY, 
strengths to zero, without changing other parameters. Then, we have recalculated the amplitudes and observables. 
This model will be hereafter referred to as "OSA, rjY excluded". The results are given in Table ITT1 and Fig. ^ The 
only observables which are not affected are the threshold ratio R c and the K~p — > 7r +v J~ cross section. All other 
quantities are significantly modified, especially 7 and R n as well as the K~p — ► K~p and K~p — > 7r~£ + cross sections 
which become unrealistic. Also, the maximum of the ttY, mass spectrum is shifted towards the higher values of the 
momentum, thus incompatible with the data. Similar effects have been pointed out in Ref. jj. This situation can 
be understood by examining the values of the strength parameters in Table IIIII The KN-rjY strengths deviate from 
unity by about 15% to 25%, and it is not possible to obtain a correct overall fit if they are constrained to stay closer 
to unity (in particular, the position of the A(1405) resonance is not correctly reproduced). This does indicate that the 
KN-rjY strengths take part in the minimization procedure at about the same level of importance as the other ones, 
and it is meaningless to turn them off entirely without re-adjusting other parameters. In order to see this problem 
from a somewhat different angle we have introduced yet another model: OSB, a variant of OSA in which the rjY 
channels arc excluded from the beginning. The values of the parameters are given in Table ITTTI The results obtained 
for the K~p branching ratios and scattering length are close to those obtained with OSA, see Table [H] This is also 
the case for the total cross sections. However, the position of the A(1405) resonance is now shifted towards lower 
values of the momentum, see Fig. 0, which means that OSB is not able to reproduce the properties of the A(1405). 
We thus need to retain the rjY channels. 

Lastly, we use the parameters thus obtained both for OS1 and OSA to calculate the amplitudes (or i-matrices) for 
the K~n initiated processes, see Eq. (|HJ|, for which, as mentioned earlier, there is no data to be confronted with. We 
still need those amplitudes for our K~d three-body calculations. Here along with the corresponding quantity in the 
K~p initiated channels, we only present the scattering lengths as found in Table llVl 

TABLE IV: KN scattering lengths (in fm) calculated in the particle basis with models OS1 and OSA. Lhe values in the last 
column have been evaluated by Ramos |60| at the same energy. a p , a° and a ex (calculated at W = M K - + M p ) are the 
scattering lengths for elastic K~p,K n, and charge exchange K~p <-> K n, respectively. a n (calculated at W = M K - + M n ) 
is the scattering length for elastic K~n. 



Reactions 


OSA 


OS1 


Oset-Ramos 


a p (K~p — > K~p) 


-0.888 + i 0.867 


-1.035 + i 0.828 


-1.013 + i 0.947 


a„(K~n — > K~n) 


0.544 + i 0.644 


0.573 + i 0.452 


0.540 + i 0.531 


aUK°n ->K°n) 


-0.444 + i 0.998 


-0.602 + i0.894 


-0.516 + i 1.053 


a ex (K~p ^K°n) 


-1.215 + j0.393 


-1.365 + i0.484 


-1.289 + i 0.484 



The differences between the results given by two models do not exceed ~ 30%. For comparison we give also the 
values obtained by Oset and Ramos p. Our OS1 results agree with those values within ~ 15%, in line with the 
differences previously observed on the other observables. The symmetry breaking effect in the mass of the hadron 
isospin multiplets is quite visible, especially in the real parts, since in the limit of isospin symmetry, one has: a p = a° , 



and a ex — a p — a n . Note that the scattering lengths in Table Hvl have been obtained at the K~p threshold (except for 
the elastic K~n process). In fact, these quantities are very sensitive to the value of the threshold at which they are 
calculated, which is then reflected in the values obtained for the A K - d scattering length. These aspects have been 
discussed in our previous paper and will be revisited at beginning of Sec. IIVI 

In conclusion, the OS1 and OSA parametrizations are good candidates to be used in the three-body calculation. 
It is clear that such interactions, the parameters of which are determined by a fit to the available observables in 
the particle basis with the chiral SU(3)-symmetry constraint, must definitely be preferred to previous separable 
interactions with parameters determined in the isospin basis, without any symmetry constraint, as it was the case in 

Refs. ESQEllIIII. 

B. NN INTERACTION 

We have considered three different relativistic separable potentials to describe the deuteron (d) channel. The 
structure of the equations may be read off from the coupled-channels cases described in Appendix El by replacing the 
particle labels with the angular momentum labels relative to the 3 <Si and 3 D\ coupled partial waves. 

All the interactions considered are of rank-1. Characteristic to such potentials, the static parameters are correctly 
reproduced, namely: the triplet effective range parameters a t and r t , the .D-state percentage value Pjj, the quadrupole 
moment Q, and the asymptotic ratio r\ = Ad /As- On the other hand, the 3 <Si and 3 Di phase shifts cannot be 
reproduced simultaneously. 

The first model (hereafter called model A) is the one proposed in our work on the nd system |61|. namely the 
parametrization denoted by SF(6.7) therein, with Pjj = 6.7%. This model is an extension of the usual Yamaguchi- 
type interactions, using form factors which are expressed as ratios of polynomials. The parameters are fitted to the 
static properties, the 3 <Si phase shift, and also to the deuteron monopole charge form factor up to about 6 fm -1 . All 
details can be found in Ref. [6l| . 

Besides this model fitted to on-shell properties only, we have considered a second model derived from the PEST1 
potential constructed by Haidenbauer and Plessas in Ref. 62]. The authors have constructed a separable represen- 
tation of the Paris potential to reproduce both its on-shell and off-shell characteristics. Among various JVJV partial 
waves, special care is devoted to the coupled 3 Si- 3 Di state. The best approximation to the Paris potential requires a 
rank-4 interaction. However, for applications where only the deuteron bound-state enters, a rank-1 parametrization 
was proposed, called PEST1, with all deuteron properties (including the wave functions) being practically the same 
as those given by the Paris potential. The price to pay is that the form factors are chosen as sum of rational functions, 
with many parameters, see Ref. [62l |. As the PEST1 parametrization is non-relativistic, we have extended it by means 
of the minimal relativity rule. More specifically, the relativistic potential V R between nucleons 1 and 2 is obtained 
from the non-relativistic one V NR according to the following transformation in momentum space: 

V R (p,p') = {2n) 3 ^j2e lp 2e 2p V NR {p, p')y/2e lp , 2e 2p , , (10) 

with €i P = y pf + to? the total energy of nucleon i. Taking V NR as separable: V NR (p,p') — Xg(p)g(p'), V R (p,p') has 
the same form, with g(p) multiplied by: 2^J (27r) 3 -\/p 2 + to 2 . We have checked that this transformation induces only 
slight changes in the deuteron properties. For example, the original Pd value is 5.8%, and the value obtained after 
the relativistic transformation is 6.1%. So, we can keep the original parameters. In our present work this is denoted 
as model B. 

Finally, in order to assess the importance of the D-state contribution in the low energy K~ d scattering, we have 
used a pure 3 S\ relativistic potential, with form factor as given in Eq. The values of the strength and range 

parameters fitted to Ejj and at are: = —5974.2, (3d — 1.412 fm -1 , respectively. This model is called model C. 

C. OTHER TWO-BODY INPUT 

We have also considered the contributions of the ttN-P^ and coupled YN interactions. As will be explained in 
Sec. IIV CI these two-body channels start to contribute from the second order in the multiple scattering expansion of 
the Faddeev equations, so their contribution in the low energy K~ d observables is expected to be moderate or even 
small. 



1. nN-Pss interaction 



We have chosen a conventional A-isobar model where the A resonance is parametrized according to a one-term 
separable potential: 

V(p,p';a)=g{ P )^-^g(p'), (11) 
a — rriQ 

Here, a is the nN total cm. energy squared, and uiq the bare A mass. A p-wave monopole form factor is assumed 
with cut-off A: 



•fp) = ^ ("> 

The strength parameter Aa, together with mo and A, are fitted to the P33 phase-shift. In the relativistic approach, 
we have obtained the following values: A A = 218.582, m = 1308.8 MeV, and A = 290.9 MeV/c. 

In order to solve the three-body equations in the particle basis, it is necessary to introduce different charge states 
of the pion and nucleon. So, the following states will contribute: the {n°p, ir + n) and {tt~p, 7r°n) coupled states, and 
the (7r - n) state. The corresponding separable two-body potentials are obtained by choosing the same form factors 
as Eq. I|12l) in all channels, and expressing the strength parameters in terms of Aa, see Appendix 1X1 The resulting 
two-body i-matrices entering the three-body equations are obtained straightforwardly. 



2. YN interactions 



The hyperon-nucleon interaction has been well studied in the past. One of the most popular approaches is meson- 
exchange potential models with 577(3) symmetry constraints used in the coupled channels equations, see Ref. jg^l 
and references therein. Besides, the intrinsic interest in investigating the available YN experimental data, these 
interactions serve as input to hypernuclear physics. Concurrently, separable models have been developed, and some 
of them have been used as input to K~d three-body calculations. In particular, the effect of the final state YN 
interaction (limited to s-waves) on the Ap invariant mass distribution near the EJV threshold was studied in Refs. 
[Til Il5| . One of the main conclusion was that the best reproduction of the shoulder in the Ap mass spectrum favoured 
models which do not support an unstable SiV bound state. The interactions that we have elaborated in Ref. 1C, 17] 
and the one presented in this work fulfill this condition. 

The data to which the adjustable parameters of the separable potential should be fitted are scarce and exhibit 
rather large error bars. To our knowledge, there are no new data in addition to those used in our previous work, 
namely: (i) the E + p — > £ + p, X~p — > £~p, S°n, An, and Ap — > Ap total cross sections, in the hyperon momentum 
range p l y h < 300 MeV/c (lA the E+p ->■ E+p, E"p -> E~p, and YTp -> An differential cross sections for p^ b ~ 300 
MeV/c. See Refs. IHUlE HE El El EllZ^I 

In Refs. 0,0], we have used most of these data to determine the parameters of two models: one non-relativistic, 
and the other relativistic. The calculations were done in the isospin basis, where one must consider the HN-AN 
(I = 1/2) coupled channels and the EA (/ = 3/2) single channel. The main conclusion was that the total cross 
sections are dominated by the 3 <Si partial wave, except for E+p — > E+p which is dominated by the 1 5'o wave, while 
the P and D partial waves have a significant contribution only in the E~p — > An total cross section (of course, these 
higher partial waves must be taken into account if the differential cross sections are added to the data to be fitted). 

The model used in the present study is an extension of the above mentioned relativistic model to calculate the 
observables in the particle basis. Only the 3 Si partial wave contributions are included, thus the parameters are fitted 
to the total cross sections, except for E+p — > E+p. Note that neglecting the 1 So partial wave is justified in the 
three-body calculation at low energies, since the contribution from the singlet S'-wave YN interactions is excluded for 
parity considerations. We take as adjustable parameters the coupling strengths and the ranges of the form factors in 
the isospin basis. The observables are calculated in the particle basis, where the following channels contribute: the 
(E°p, E+n, Ap) and (E~p, E°n,An) coupled channels, and the E~n and E+p single channels. The relations between 
the strengths parameters in the two bases can be found in Appendix EI an d the transition matrices for the different 
reactions are obtained from the general expressions given in Appendix 151 We take monopole form factors, Eq. 
with isospin- independent ranges. The values of the fitted parameters are given in Table IVl and the cross sections 
in Fig. |21 The selected data can be well reproduced by this simple model, but a large reduction effect is observed 
when isospin basis is adopted in the total cross sections for the S _ p induced reactions at low values of the hyperon 
momentum. 




FIG. 3: E"p -> E"p, E°n, An, and Ap — > Ap total cross sections, calculated in the particle basis (full line) and in the isospin 
basis (dashed line). Experimental data are from Refs. [64L Ififll. Ififl. l67t l68t IrM l7flj . 



TABLE V: Strength and range parameters of the YN interactions. 



Isospin Channel Parameters 

1/2 EiV -3209.69 

ZN-AN -1739.22 

AN -1794.16 

3/2 EiV 3072.26 

EiV range (MeV) 452.300 

AiV range (MeV) 356.981 



We give in Table Wl\ the YN scattering lengths calculated in the particle basis. 

Except for a^+p-s+p and o,s- n —s- n which are practically equal, the symmetry breaking effects are large. For exam- 
ple, at the limit of exact isospin symmetry, one should have: as-p-s°n = V^{o.s°n-'E o n — a £-p-£-p) — ~ a s+n-£p = 
°p-s°p — a s+n-s+n)j which is clearly not the case. 



TABLE VI: YW scattering lengths (in fm) calculated at W = M E+ + M n for channels (E + n, E°p, Ap), and W = M E - + M p 
for channels (E~p, E°n, An). The single channels values are : a s - n = —0.454 (W = M E - + M n ), and a s + p — —0.455 
(W = M E+ +M P ). 
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1.851 



III. THREE-BODY EQUATIONS FOR THE K~d SYSTEM 



In this Section, we describe the three-body equations for the K~d system, in which the two-body input described in 
the previous Section will enter. In the first Subsection, we recall the general form of the three-body equations written 
in the isospin basis, in the case of coupled two-body input channels, and we give the equations for the rotationally 
invariant amplitudes. Then, the antisymmetrization due to the identity of the two mucleons in the isospin basis is 
examined. In the second Subsection, the extension to the particle basis is given. Some important aspects concerning 
the practical calculation are considered in the third Subsection. 

A. Three-body equations in the isospin basis 

The extension of the usual three-body equations to the case where the two-body operators connect two states 
involving particles which are different (i.e. inelastic channels), results in the following system of coupled equations, 
written in operator form: 

X ab (s) = Z ab (s) + Z ac { S )R cc ,{s) X c , b (s), (13) 

c, c! 

with s the 3-body total energy. Here, a, 6, c and d are the indices which specify the particles involved in the various 
three-body channels, namely the spectator and the interacting pair. X ab is the transition amplitude between channels 
a and b, and Z ab is the corresponding Born term. The main difference with respect to the usual case is that the 
two-body operator R cc i connects two different two-body states labeled as c and d . 

We now specify the values taken by the channel indices in Eq. (|13|) . Taking into account all the two-body input 
considered in Sectional one must consider the following types of three-body channels in the isospin basis: d [K(NN)], 
y[N{KN)}, a[N(wY)], fi[N(r)Y)], (3[n(YN)}, and A[Y(7riV)]. Here, the first letter is the label of the channel, 
and in square braces we specify the spectator and the associated pair of particles. In fact, labels y, a, [i and f3 can 
be considered as "generic" names. In practice, extra indices are needed to fully describe the physical situation. For 
example, for the pairs corresponding to the KN interactions, there are three 1 = coupled channels (KN-ttT,-t]A), and 
four 1=1 coupled channels (KN-nTi-TrA-riA)). A similar situation occurs in the case of the (3 channels corresponding 
to the T,N-AN (I = 1/2) coupled channels. The corresponding three-body channels are labeled as summarized in 
Table IVTT1 Note that in channel A, the spectator hyperon Y is restricted to S, from isospin considerations. 



TABLE VII: Labels of the three-body channels. The second line specifies the isospin of the two-body sub-system. 



channel K(NN) 


N(KN) 


A(ttS) 




N(KN) 


A(ttE) A(ttA) 




S(ttA) 


tt(EA) 


ty(AN) 


tt(SA)) 


isospin 











1 


1 1 


1 


3/2 


1/2 


1/2 


3/2 


label d 


yi 


Ol 


Mi 


J/2 


«2 as 


M2 


A 


ft 


ft 


ft 



In the formal equations (|13f) . the channel indices a, 6, c, d take their values in the set defined above: {d, y, a, /i, f3, A}, 
and the different quantities: X, Z, R, can be considered as matrices with respect to these indices. 

The transitions between different two-particle states are induced by the two-body operators R CC ''- for example, 
R ya is the quantity R^n^^y corresponding to the two-body t-matrix for the KN — > nY transition, evaluated in the 
presence of a spectator nucleon. Written in matrix form, the non-zero R operators appear as block-matrices as shown 
in Table IVHll 

Concerning the Born terms, only those which connect initial and final states involving the same three particles are 
different from zero. These terms are shown in matrix form in Table llXl 

For example, Zd y is the Born term for the exchange of a nucleon between the K(NN) and N(KN) states. Note 
that Zdd is equal to zero, since no particle can be exchanged between the initial and final NN pairs. We note also 
that, using the two-body input channels considered here, there is no Born term involving an r/Y pair. To have such 
terms, it would be neccessary to take into account the contributions of three-body channels like rj(YN), Y(wN) and 
Y{qN). (the last two channels necessitate to introduce the ttN-tjN two-body coupled system). These contributions 
are expected to be negligible at the K~d threshold. 

The successive steps leading from the formal equations Ijl3|l to the relativistic equations for the rotationally invariant 
amplitudes are the same as in the usual case. We refer the reader to References 0,0>EHEllz2 for all details. The 
final equations read: 



TABLE VIII: Matrix of propagators in the isospin basis. 
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TABLE IX: 


Matrix of Born terms 


in the isospin 


basis. 


For the 


non-zero 


Born terms, the exchanged particle 


is shown in 



parentheses. 
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X^> ,p c ;s) = 3£ c (p ,p c ;s) + J2 J ^Zil(p a , Pb - S )Rll7 Cb '(a b )X^(P^Pc; S ), (14) 

b,T b ;b',T b i 

where a b is the invariant energy of the pair in channel b expressed in the three-body center of mass system, c a = 
(J a ,S a ,I a ) specifies the conserved quantum numbers of the pair in channel a, and r a = (c a ,^ a ,E a ) specifies the 
three-body quantum numbers in this channel. Note that labels c and r refer to the spin-isopin variables in a given 
channel. For example, assuming that channel a is composed with particle i as spectator and the pair (jk), we define 
the following quantities: 

- Si : spin of particle i , 

- Si (= Sj + Sfc), Li, and Ji (= Li + Si) : spin, orbital angular momentum, and total angular momentum, 
respectively, of pair (jk) , 

- Si (= Si + Ji), li, and J (= h + Xj) : channel spin, orbital angular momentum of i and (jk), and three-body 
total angular momentum, respectively. 

The isospin variables are defined in the same way: 

- ti : isospin of particle i , 



- Ii (= tj + tk) ■ isospin of pair (jk), 
X (= ti + Ii) : three-body total isospin. 

The notations used in Eq. (|14|l for the two-body propagators are the following: the lower indices b and b' refer 
to the involved coupled channels, and the upper "index" Cf, = q,< means that the spin-isospin quantum numbers are 
conserved by the interaction, i.e.: 3\, = Jy, Sb = Sv, h = h' (the latter equality holds only in the isospin basis). For 
the uncoupled propagators, the single index C{, specifies completely the interacting pair. 

The two-body propagators are calculated as explained in Section Inland Appendix iBl and the general expression of 
the Born term can be found in Refs. |l7j . 

We end this Subsection with two specific aspects concerning the calculation of the Born terms in the isospin basis. 
The first one concerns the problem of antisymmetrization. In the isospin basis, the particles in the different multiplets 
are considered as identical. In particular, the neutron and proton are treated as identical particles: the nucleon N. 
Therefore, one must properly antisymmetrize the amplitudes and Born terms where the initial and final three-body 
states involve two nucleons. As a result, one must introduce antisymmetrization coefficients as explained in Appendix 
ICl The second aspect concerns the problem relative to the "ordering" of particles when evaluating the Born terms. 
For example, let us consider the Z yy Born term for K exchange between two KN states, with a spectator nucleon. 
The antisymmetrized expression is, according to Appendix [U] 

Z yy = -Zyty* = -(N 2 (KNi) \G \ Ni(KN 2 )), (15) 

with Go the three-body propagator. Usually, the expressions given for the Born terms assume a cyclic ordering of 
the particle labels, see for example Refs. 0,|zil Namely, one calculates for example: = (i(jk) \Go\j(ki)), for 
the exchange of particle k between pairs (jk) and (hi), where i, j, k are assumed to be cyclically ordered both in the 
initial and final states. In the case of Z y i y 2, Eq. (|15(l . the spectator particles in the initial and final states are labeled 
as: N 2 = i, N\ = j, thus the Born term has the "non-cyclic" form: (i(kj) \Gq\ j(ki)) . We obtain the cyclic form 
by exchanging particles j and k (i.e. N and if), in the pair of the final state. This introduces the following phase 
coefficient: 

(_!)%« (_i) s i?iv- s «- ;5 ir(_i)%iv- f «- t i? . (16) 

The first factor is due to changing the direction of the relative momentum of the pair, and the second (third) factor 
results from the property of the Clebsch-Gordan coefficients when the coupling order of the spins (isospins) of two 
particles is changed. 

B. Extension to particle basis 

In the particle basis, we consider the deuteron as composed of two distinct particles: the neutron (n) and the proton 
(p), and the particles of the different multiplets take their physical masses. The number of three-body channels to be 
considered increases considerably as compared to the isospin basis case. For example, we have seen in Section III Al 
that for the KN interactions in the particle basis, we must consider the coupled channels related to K~p, namely: 
K~p,K°n, tt-Y,+ , 7T+IT, 7r E , vr°A, 77s , 77A, and those related to R-n: K~n, tt^E , 7r°Xr , 7r _ A, ^E". So, if wc 
retain only these contributions in addition to the deuteron, the three-body channels to be considered are: the K~ (pn) 
channel, the eight above channels with the neutron as spectator, and the five remaining channels with the proton as 
spectator. So, we have a 14x 14 Born terms matrix. However, the number of non-zero Born terms is very limited. In 
the case considered here, we have only the following different Born terms: ( K~(pn) \Gq \ n(K~p) ) (proton exchange 
between the deuteron and the (K~p) pair), (K~(pn) \Go\p(K~n) ) (neutron exchange between the deuteron and the 
(K~n) pair), ( n(K ~p) \Gq\ p(K~n) ) (K~~ exchange between the (K~p) and (K~n) pairs), (n(K n)\Go\n(K n)) 
(K exchange between the initial and final (K n) pairs), and the symmetric terms. 

The matrix of propagators has a block structure similar to that in the isospin basis (Table IVIIII ) : single term for 
the deuteron propagator, 8x8 matrix for the channels coupled to K~p, and 5x5 for the channels coupled to K~n. 

If we take into account the contributions of the ttN and YN interactions, the following additionnal three-body 
channels must be considered (see Section III C|) : {E _ (7r°p), £~(7r + n)}, {E° (ir~ p) , Y,° (ir° n)} , and S + (7r~n), coming 
from the nN interaction, and: {tt~ (S + n), it~ (£°p), ir~ (Ap)}, {tt (£,~p), 7r (E n), 7r°(An)}, and 7r + (S _ n), coming 
from the YN interactions. 

The derivation of the equations for the rotationally invariant amplitudes, and the calculation of the Born terms, 
are done along the same lines as in the isospin basis. The main difference is that there is no problem relative 
to antisymmetrization, except for (n(K n)\Go\n(K n)) which must be antisymmetrized with respect to the two 



identical neutrons, as in the case of the Z yy Born term calculated in the isospin basis, see Appendix El Another 
important difference concerns the isospin dependence. In the isospin basis, the expression of the Born term involves a 
"6-j" coefficient originating from the transformation from initial to final three-particle isospin states, see Ref. 
This coefficient depends on the values of the particle isospins, total isospin of the initial and final pairs, and three-body 
total isospin. In the particle basis, the individual isospins of all particles are well defined, but not the total isospin of 
the pairs (consequently, the c labels in Eq. Ijl4|l do not depend on the isospin I of the pairs). As the initial and final 
three-body states involve the same three particles, the isospin coefficient simply reduces to unity. 

C. Practical calculation 

For the practical calculation, we must first define the values of the various isospins, spins and angular momenta. At 
first, we consider the isospin basis case. The total isospin of the K~d system is I = 1/2. Now, in a given three-body 
channel, the quantum numbers (L, S, J, 7) of the pair and the spin s of the spectator particle are fixed. The channel 
spin S is then given by: \s — J\ < £ < s + J, and the angular momentum I of the spectator relative to the pair by: 
\J — E| < I < J + E, with J the three-body total angular momentum. For a given value of J7", the possible values of I 
can be ordered in two sets corresponding to opposite parities of the three-body system. The situation is summarized 
in Table K| 

TABLE X: Two-body (L, S, J, I) and three-body (I, E, J) quantum numbers in the isospin basis. The two-body partial waves 
are labeled as : 2S+1 Lj for NN and YN, Z/2/.2J for nN, and Li^j iorKN and txY . The column labeled as l a (h) corresponds 
to negative (positive) parity states for odd values of J, and to positive (negative) parity states for even values of J. Only the 
values I > are retained. 
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In the present paper, we consider only the K d scattering length which is defined as: 

A K-d = - lim \ r X dd , (17) 

where Xdd is the (J = 1~, / = /' = 0) partial amplitude for K~d elastic scattering, evaluated at the zero limit for 
the kaon momentum. If we retain the contributions of the d+KN+A+YN two-body channels, we have a system of 
12 coupled three-body channels (see Table IvTTjl . After angular momentum reduction, we obtain for J = 1~ a system 
of 25 coupled equations, see Table Ixl The singularities of the kernel ZX are avoided by using the rotated contour 



mpic 

method [lOj . and, after discretization of the integrals, this system is transformed into a system of linear equations. 



It is well known that the iterated form of the three-body equations does not converge. This is illustrated with the 
following results obtained (in the isopsin basis) with the "OSA+deuteron-A" model. The values of A K - d obtained at 
first, second and third order of iteration are respectively (in fm): (0.303 + i 3.258), (—0.572 + i 3.684), and (—1.366 + 
i 4.131), which clearly do not converge. To obtain the exact value, we solve the linear system by matrix inversion. In 
the above example, we get: (—1.636 + i 2.618) fm. We can also use the Pade approximants technique which leads to a 
convergent solution from the successive iterated terms. In practice, we have used a diagonal [5/5] Pade (constructed 
with the 11 first iterates), which was found to be sufficient to achieve convergence. The dimension of the matrix to 
be inverted is rather large, especially when the contributions of the "small" two-body partial waves are taken into 
account. However, this is a sparse matrix because of the limited number of non-zero Born terms. Thus, it is much 
less time consuming to solve the linear system with using the Pade approximants method. 

The extension to the particle basis case is straightforward. As explained above, in this case the Pade approximant 
method must be even more prefcred to the usual methods for solving the linear system, since the number of coupled 
channels is much larger than in the isospin basis for a given choice of two-body contributions, but with a large number 
of Born terms equal to zero. 

Note also that we have checked that, using the particle basis computer program with the particle masses in the 
multiplets replaced with the mean values used in the isospin basis, we have obtained again the isospin basis results. 



In this Section, we present our result for the scattering length A K - d . As an introduction, we argu e why we think it 
necessary to go beyond the Fixed Centre Approximation (FC A) as adopted by Kamalov et al. [41| . Then we discuss 
some general aspects regarding the choice of the basis (viz. particle vs. isospin) and the use of different coupled 
two-body K N input, and we investigate the effects due to the choice of different deuteron models. At the end, the 
effect of "small" two-body input on the K~d scattering length is discussed. 



In the Introduction, it has been stated that we need to go beyond the FCA and to solve the three-body equation 
exactly. We justify that claim here. 



IV. RESULTS 



A. 



Why need to go beyond FCA ? 
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FIG. 4: K p elastic scattering amplitude as a function of the cm. total energy, calculated in the particle basis with the 
"OSA+deuteron-A" model. Full line: real part, dashed line: imaginary part. 



In the FCA the deuteron is viewed as composed of a proton and a neutron with a fixed separation r. The incoming 
zero energy if "-meson then makes multiple scattering off the proton and neutron with no recoil of the target particles. 
Within this approximation the three-body scattering equation may be solved algebraically: see Eq. (23) of ref. |4lj . 
to find the scattering length operator A K ~ d (r) expressed in terms of the KN scattering lengths: a p , a n , a° , a ex 
as in Table llVl and the separation r. The actual K~d scattering length is then identified as the expectation value 
< ipd\Aic- d,\4>d > over r with respect to the deuteron wave function ip d (r). So essentially, A K - d is determined by the 
two-body if A" scattering lengths mentioned above. 

As discussed in Section |n] the A(1405) (7 = 0) resonance is generated as a bound state of K~p embedded in 
the irY continuum. Now the position of this resonance is fairely close to the threshold for K~p (w 1432 MeV) and 
if n (ss 1437 MeV), respectively. Thus the elastic K~p, charge exchange: K~p — > if n, and hence the elastic if n 
scattering are all affected by this resonance, and the corresponding amplitudes vary rapidly near their thresholds: 
see Fig. 0] for the case of the elastic K~p amplitude which shows strong variations, particularly in its real part. 
The question then arises as to at which energy these amplitudes should be calculated to produce the corresponding 
scattering length in use for calculating A K - d in FCA. The binding energy of the deuteron is w 2.25 MeV. In Fig. 03 




FIG. 5: KN and K d scattering lengths calculated in the particle basis for the following values of the two-body threshold 
energy W: M K - +M P = 1432 MeV, M K - +M n = 1432.65 MeV, M K - +{M p + M n )/2 = 1433.3 MeV, and A%o + M n = 1437.3 
MeV. Model OSA is used for the KN interaction. The K~d scattering length is calculated with the FCA approximation, using 
the " OSA+deuteron-A" model. Symbols are placed at the threshold values. The lines are to guide the eyes. 

we show the values of the two-body scattering lengths which enter the FCA calculation and the resultant A K - d for 
different values adopted for the KN threshold W . We immediately notice that because of the strong variation of 
the two-body scattering lengths, the corresponding A K - d varies also rapidly with a slight shift in W: the real part 
in particular. This demonstrates that just because of the proximity to the A(1405) resonance, the FCA is not very 
reliable. In a full three-body results, the quantity corresponding to W is ^fa: the energy available to any two-body 
amplitude in the presence of a spectator. This is smeared out due to a loop momentum integration. Consequently, 
the full three-body results do not suffer from this excessive sensitivity. 



B. Calculations including only the if TV and NN interactions 

Here we retain only the deuteron and the if A two-body input in the three-body calculation of the K~d scattering 
length. This means that in the multiple scattering expansion of the equations only the following three-body states 
do enter: K~{pn), n(K~p), n(K n) and p(K~n). Although no explicit three-body states with hyperons enter in 



the calculation in this approximation, there is an implicit effect from the irY channels through the two-body KN 
input. Hence it is more than the single channel approximation studied in Ref . |74| . Within the present context we 
first confront the results from the particle and isospin bases, which makes us keep the result with particle basis in the 
subsequent subsections. We then study the results with different KN models. 

1. Isopin and particle bases 

The isospin symmetry breaking effects in the KN sector have been clearly demonstrated in Sec. Ill Al by comparing 
the K~p observables obtained both in the isospin and particle bases (see Table ITT1 and Figs. ^and[3J). In particular, 
upon going from the particle to isospin basis, we find that the magnitude of the real part of a^~ p obtained from 
model "OSA+deuteron-A" decreases by about 20% whereas the imaginary part increases by as much as 30% (see 
Table ITT)l . This last tendency should be due to the fact that in the isospin basis (we reiterate here that in this case 
the isospin symmetry is exact) two-body channel thresholds become identical among different KN channels within 
the same iso-multiplets. Consequently, these effects are reflected in A K - d , as one can see Table fxTI for different KN 
models: by going from the particle to isospin basis, the magnitude of its real part gets smaller by about 10%, while 
its imaginary part gets strongly enhanced by ~ 70%. In what follows we will mostly comment on the results from the 
particle basis since they are more realistic. 

2. Results using different KN models 

We give in Table IXll the values of A K - d obtained from models OS1, OSA and OSB of the K N interactions described 
in Section Til Al We briefly recall that by construction all the models have adopted the SU(3) symmetry constraint 
on the strength parameters, with the possible breaking effects introduced in OSA and OSB, and the contributions of 
the rjY channels have been taken into account, except for OSB. 



TABLE XI: Sensitivity of the K d scattering length to the model used for the KN interaction. The deuteron is model A. The 
calculations are done in the particle and isospin bases. 



Model 


OS1 


OSA 


OSB 


particle basis 


-1.985 + il.642 


-1.802 + i 1.546 


-1.722 + i 1.354 


isospin basis 


-1.759 +i 2.907 


-1.636 + i 2.618 


-1.709 + i 2.247 



The calculations have been done with model A for the deuteron. The variation observed on the real and imaginary 
parts of A K - d calculated in the particle basis are moderate: for example, the real part increases by about 9% and 
the imaginary part decreases by about 6%, when replacing model OS1 by model OSA. The absolute values of these 
variations are comparable to those observed in the K~p scattering length (see Table ITT]) : by going from OS1 to OSA 
we find an increases in the real part by about 13%. For the imaginary part the corresponding increase is less than 
1%. 

Concerning the contribution from the rjY channels in the KN interactions, we observe the same type of effects, 
both on a^-p and A K - d . Specifically, compared with the values given by model OSA, both the real and imaginary 
parts of a-K-p decrease by ~ 9% when using OSB (see Table ITTj) . while the real part of A^-d increases by ~ 4% and 
the imaginary part decreases by ~ 12%. Now we recall the result in Section Til Al using model OSA and excluding 
the contributions of the rjY channels, we have seen that all KN observables were strongly affected. In particular, the 
imaginary part of a K - p increases by as much as ~ 60%, and the real part by ~ 18% (see Table ITljl. A similar tendency 
is observed on A^-d, as the OSA value of (—1.802 + i 1.546) fm has been transformed to (—1.705 + i 2.360) fm when 
the r\Y channel contributions are excluded (as discussed in Section HTAl this means that the coupled channels equations 
were solved again with the parameters of OSA but excluding the coupling to the rjY channels). In conjunction with 
the observation in Section fll Al iust mentioned above, we conclude here that it is important to retain these rjY channels 
for Apc~d- Before ending this subsection we should mention that, as shown in our previous publication |l|, we have 
also tested the KN amplitude from Ref.[j| in the three-body equation, and the result turned out to be very close to 
the one with OS1, as expected. 



3. Dependence on deuteron models 



We have tested three different deuteron models: A, B, C as described in Sec. Ill Bl Their D state probabilities: Pjj, 
are 6.7%, 5.8%, and 0%, respectively. These deuteron models are used in combination with the OS1, OSA and OSB 
parametrizations of the KN interaction, and results are summarized in Table IxTTl 



TABLE XII: Sensitivity of the K~d scattering length, calculated in the particle basis, to the model used for the deuteron 
channel. Models A and B have D-state percentage values of 6.7% and 5.8%, respectively, and model C is pure 3 5*1. Models 
OS1, OSA and OSB of the KN interactions are considered. 



Model 


A 




B 




c 




OS1 


-1.985 + ' 


il.642 


-1.966 +-, 


•; 1.515 


-1.975 + 


i 1.313 


OSA 


-1.802 + 


i 1.546 


-1.788 + -, 


i 1.435 


-1.780 + 


i 1.243 


OSB 


-1.722 + : 


i 1.354 


-1.703+-, 


i 1.263 


-1.685 + 


i 1.128 



Irrespective of which KN model is adopted, there is a definite pattern: the imaginary part of the scattering length 
increases in accordance with the corresponding increase in the _D-state probability of the deuteron. The increase in 
this quantity due to a change: Pn = — > 6.7%, is as large as 20%. The real part appears to stay more or less the 
same in the meantime. We have looked at a few first terms in the three-body multiple scattering series, but that does 
not tell us why the change is almost exclusively in the imaginary part. But whether it is in real or imaginary part, 
the following simple observation should suffice in understanding the change of this magnitude. First, we introduce 
the S and D component of the deuteron wave function in momentum space as tps(p), i^D^p), such that 

/>oo />oo 

P s = ^ 2 s (p)p 2 dp, P D = i> 2 D (p) P 2 dp, P S + P D = 1, (18) 
Jo Jo 

where, for example, we may take Ps ~ 0.93, Pu ~ 0.07 for model A deuteron. We may then re-write 

i>s($) = \[P~s4>s{v\ ^d{p) = Vp^Mp), (19) 

such that both <frs and <f>n are normalized to unity. Note here that \fP~s ~ 0.97 and vPd ~ 0.26, respectively. Then 
we may express the scattering length in the following manner: 

A K -d = Ps< <Ps\A SS \<Ps > +2y / P D P s < <t>D\A D s\<t>s > +Pd < tf>D\A DD \fo >, (20) 

where A are operators in the space of the deuteron wave function, and the closed brackets mean the integration over 
the initial and final loop (or off-shell) momenta. In the above expression, it is easy to see that due to the second S-D 
interference term the result with and without the deuteron D state could differ up to a few 10's of a percent, even 
though Pjj is just about ~ 6%. Of course it is likely that the matrix element < <j>D\ADs\4>s > m ay well be smaller 
than the one for the first term: < 4's\Ass\ ( l : 'S >; due to a slight angular momentum mismatch between the S and 
D states. However, when we study the behavior of ips(p) and ipn(p), it is easy to observe the following trend: the 
former is very large at p — and decreases rather rapidly down to p = 0.6 fm _1 . The latter is zero at p = to 
start with, but it's magnitude increases rapidly up to about p = 0.6 fm _1 , then decreases moderately down to about 
p = 2.5 fm _1 . The result is that the wave function components become of the same order of magnitude from about 
Pequai = 0.75 fm _ . On the other hand, we find that p eqU ai is in the sub-threshold region for the K~ p amplitude 
dominated by the A(1405) resonance, see Fig. 0] This is imbedded in Ar>s- For this reason it may be fairly likely that 
< 4>d\Ads\4>s > and < 'pslAssl'Ps > are not very different in magnitude, hence a ~ 20% increase in the imaginary 
part of A K - d due to the deuteron D state is possible. So it is important to retain that component. 

C. Effects of the "small" two-body input 

What we term here as small two-body input are (i) the irY channels resulting from the initial KN interaction, (ii) 
ttN appearing with a spectator hyperon, and (iii) the YN interactions with a spectator nucleon in the three-body 
equations. Here we exclude the channels involving the r\ meson, recall Tables IVIII In this manner, in principle, our 



equations do satisfy three-body unitarity exactly. It was shown in Refs. |l4l Il5| that these interactions were very 
important in the threshold break-up reactions: K~d — > TtNY as they control the final state interactions. Here we 
want to study these effects in the threshold elastic case. For some convenience, calculations in the isospin basis come 
back in our discussion. 

We take models OS1 and OSA for the KN interactions, and model A for the deuteron. Then, we add successively 
the •kN-Pzz and YN interactions described in Sec. Ill CI The results are given in Table IXTTTl the third and fourth 
columns give the values obtained when only the ttN-P^ or YN input are added to the KN+d input, and in the last 
column both contributions are taken into account. 



TABLE XIII: Contributions of the small two-body input to the K d scattering length (values in fm). Models OS1 and OSA 
are used for the KN interaction, and model A for the deuteron. The calculations are done in the particle basis. 













Model 


d+KN 


+ A 


+ YN 


+ A + YN 


OS1 


-1.985 +i 1.642 


-1.985 + il.663 


-1.975 + i 1.611 


-1.974 + i 1.634 


OSA 


-1.802 +i 1.546 


-1.793 + i 1.562 


-1.805 + i 1.511 


-1.796 + i 1.529 



The effect of the additionnal two-body input is negligible. This can be easily understood by explicitcly writing the 
three-body equations. For example, let us consider Eq. I|13(l written in the isospin basis, with the following simplified 
two-body input: d +KN (limited to / = 0, without the rjA channel contribution) + ttN-P^. Using the channel labels 
as defined in Table IVTT1 (with y\, ot\ simplified into y, a, respectively), the explicit form of the coupled three-body 
equations is: 

Z dyRyyXyd T Z dyRyaX ad 

Zyd T ZydRdXdd ~t~ ZyyRyyXyd T ZyyRy a X a d (21) 
ZaARA^Ad 

ZAaRayXyd + Z AaRaaX a d 

Using the last equation to express X a d in the first three equations, we see that the ttN-P^ channel contributes in 
terms of second or higher order, thus its effect on the low energy K ~d observables like A K - d should be small, although 
the resulting ttN-P^ state in the absence of a spectator Y is nearly on the A resonance peak in its two-body center 
of mass energy. This is consistent with a semi-quantitative estimate of the effect by Kamalov et al. Similar 
arguments hold for the contribution of the YN interactions, and also when the particle basis is used. 

Now, we may need to discuss the problem associated with the fact that the signs of the off-diagonal parameters of 
the KN and YN interactions are undetermined. Let us consider for example the OSA model of the KN interaction. 
As explained in Section Ul Al the signs of the off-diagonal strengths are those of the SU(3) coefficients given in Tables II 
and III of Ref . [f| . Now, if these signs are changed, the signs of the corresponding off-diagonal two-body propagators 
are also changed, but the KN observables are not (except the signs of the corresponding scattering lengths). For 
example, changing the signs of ^kn-ttH both for 1 = and 1 = 1, and/or the sign of Ag=jy_ 7rA (which contributes only 
for I = 1), does not affect the observables. Similar conclusions hold for the YN interactions when the sign of \sn-an 
is reversed. Now, the situation is not so simple in the three-body sector. To examine what happens, we consider the 
following cases in the isospin basis (we choose the isospin basis for the sake of simplicity in handling labels, but the 
conclusions are the same in the particle basis): 

i) Only the d+KN interactions (without the rfY channels contributions) are adopted. The corresponding three-body 
equations have the form J2J, with Z a A = Zao, = 0. Thus Eqs. I|21() reduce to the first two equations with X a d = 0. 
As these equations do not involve the off-diagonal KN-ixY propagators, the reaction amplitudes do not depend on 
the sign of the off-diagonal KN-irY coupling constants. 

ii) Next, we add the iiN-P^ two-body channel to the one just mentione above. Going back to Eqs. I|21|) . it is clear 
that we obtain the same sytem of equations if we change the signs of both the off-diagonal KN-ttY propagators and 
the X a d and XAd amplitudes. Therefore, the K~d scattering length will not be affected when changing the signs of 
the off-diagonal KN-ttY coupling constants. 

in) Now, we add YN interactions to the model with the d+KN interactions, and change the sign of the EJV- 
AN coupling constant. Then, contrary to the previous case, the K~d scattering length gets changed. This can 
be understood by noting that, due to isospin conservation, the A exchange in the three-body sector is possible only 
between the tt(AN)i =1 / 2 and N(ttA)j = i states (Born term Zp 2Cl3 , see Table llXj) . This "dissymmetry" (comparing with 
the situation for the £ exchange) implies that we cannot change simultaneously the signs of the T,N-AN propagators 



Xdd — 

Xyd = 

X a d — 

XAd — 



and of some of the three-body amplitudes without changing the original system of equations. However, as explained 
before, the change in A K ~ d does not exceed a few percent (see Table IXlTT)) . therefore this problem will not be regarded 
as an important issue worthy of extensive discussion in the present paper. It should be appropriate to stress in this 
regard that we do not anticipate any significant lack of precision because we have adopted a separable rank one form 
for the hyperon-nucleon 5- wave interactions (which are the part of the "small" input): we have clearly witnessed that 
these channels have been found to give only a small effect in the calculation of A K - d . In particular, if we accept 
that the YN interactions are dictated by SU(3) symmetry, just like om KN two-body input, and thus adjust all the 
signs of strengths in our separable potentials, for example, to the corresponding S'-wave projected Nijmegen meson 
exchange potentials [6j|, then the sign ambi guit y will also be gone out of our discussion. In fact the related problem 
was already studied by Dalitz et al., Ref. |l5j |. who found important variations in the Ap mass spectrum in the 
threshold break up reaction: K~d — > irNY. This will necessitate us to re-examine the break up channels within our 
present approach. At present, we assess that the effect of the "small" two-body input is actually not important for 
the calculation of A K - d . 

V. DISCUSSION 

In this Section we present the best estimate for the theoretical value for A K - d in our three-body approach. Then 
we discuss what kind of uncertainty should be associated with this value which derives from the effects we have left 
out in the present calculation. 

From the result presented in the last Section, we make the following choice for our preferred value: the one from the 
combination of the two-body input "OSA+deuteron-A" , in the particle basis, see Tables IXll and IXTfl for the all set of 
the results. We have not adopted the values from the calculation explicitly incorporating the hyperon channels which 
are associated with what we called small two-body input in Sec. [HJ as should be clear from Table IXIIII the effects are 
found to be quite small (< 3%). Besides, effects due to the sign ambiguity in the off-diagonal YN amplitudes do not 
exceed 1%. So we take a conservative estimate of the K~d scattering length to be 

A K - d = (-1.80 + i 1.55) fm, (22) 

to which we may assign a possible uncertainty of a few percent. 

We then need to assess the effects which may not have been taken into account in an ordinary three-body equations 
approach like the present one. The first such processes are possible four-body intermediate states: those with two 
mesons and two baryons. Diagrammatically, they may be divided into partially and totally connected ones. Of 
partially connected diagrams, those associated with baryon self energies should be dropped from consideration since 
we assume to have been dealing with the initial Kd channel. Then the remaining partially connected diagrams are (i) 
the ones in which there is a spectator meson and two baryons exchanging a meson, and (ii) the ones with a spectator 
baryon, and a meson and a baryon exchanging a meson. For the first ones they have already been included effectively 
in the input baryon-baryon (NN or YN) interactions. Likewise, the second ones are effectively included in the coupled 
KN and ttN input since they have been fitted to data. So we have only to worry about the completely connected 
diagrams. Quite fortunately, except practically for a couple of diagrams 3 , they are reduced to two baryon (YN) 
interactions: crossed two meson exchanges, and one meson exchanges with meson-baryon-baryon vertex corrections 
due to virtual meson creation and absorption across the vertex. So those completely connected diagrams are just the 
pure two-body intermediate YN channels resulting from absorbing if (or ir). 

We thus should consider only the effects of the meson (K~) absorption in the K~d elastic scattering at threshold. 

In Ref. [TtJ the p-wave effect in the low energy KN interaction has been studied within the context of an effective 
chiral Lagrangian. In that work, the effect derives from the s-channel pole contributions (the absorption/re-emission 
of if by a nucleon: KN — > Y — » KN). So this could be used as a measure for the meson absorption effect under 
consideration. The authors have stated that the p-w&ve effect is quite small: total cross sections change very little, 
whereas the differential cross sections have improved to follow the trend seen in experimental data. From this 



3 Those may be classified as meson exchange current terms that are found as diagrams (b) and (c) in the work of Wei nbe rg |7Rj for the 
case of the pion-deuteron scattering length. This combination has been calculated earlier by Robilotta and Wilkin l76l and found to 
contribute less than 3% to the total scattering length. Weinberg's calculation did confirm this earlier result. We may conclude that 
contribution from the corresponding diagrams for A K - d be far smaller. This is based upon the forthcoming discussi on in the main text, 
on the comparative study of A K - d and A n - d related to the meson absorption effect. Note that as explained in Ref. |7Sj. a combination 
of diagrams (d) and (e) should give no contribution for the scattering length. 



publication, what we could possibly exploit as the indicators of the K absorption effect semi-quantitatively may be 
the bare mass of A, viz. M\ or the ratio R c , recall Section [H] for discussion on the available experimental data for 
the coupled K~p channels. In this regard, as found in Table ITT1 the other two ratios: 7 and R n are too sensitive to 
be used for our objective. We have observed that the change in R c by the K~ absorption effect (viz. the inclusion of 
the p-wave in the language of Ref.[73) is at most ~ 2%, whereas the shift from the bare mass to the physical one due 
to the same effect for the A is about 3% (the corresponding values following Eq.(26) in that publication cannot be 
used since the basis model for the s-waveKN interaction has been modified by readjusting the subtraction constants 
in Eq. (26)). Since the kinetic energy available to the YN system after the K~ absorption in the K~d system at 
threshold corresponds to p l y b > 1 GeV/c, the effect of the YN interaction is expected to be small (a reasonable guess 
may be reached from Fig. 3 in the present article and in Ref. 63] ). In such a situation what we have just estimated 
above may well be interpreted as the effect of the K~ absorption. To be on the conservative side we set this to be a 
possible correction of a few percent. 

Not directly applicable but rather useful information regarding the effect of meson absorption comes from the 
Tt~d scattering length: A^- d . With the exception of some small effects from tt~ d — > ix°nn and n~d — > jnn, the 
scattering length in this process is purely real if no strong pion absorption effect is in effect, viz. no imaginary part 
in the absence of pion absorption. With this in mind, ealier model calculations indicated that the effect creates 
contributions both to real and imaginary parts of the scattering length. An earlier three-body model calculations 
[78j obtained TZe(A 7r - d ) = —0.035, and Im(A 7T - d ) to be between 0.0062 and 0.0075, both in units of the inverse pion 
mass: m~ x . What should be emphasized here is that the pion absorption contributes to the real part (commonly 
termed as the dispersive effect) with just about the same magnitude as the imaginary part, but with the negative 
sign. This characteristic feature was confirmed by a later calculation in multiple scattering in Ref. |79| . and the pion 
absorption contribution was evaluated to be 

AAff d = (-0.008 + i 0.011) rn,- 1 . (23) 

A more complete three-body calculation explicitly including the pion absorption in a fully consistent manner |scj 
obtained 

A w - d = (-0.047 + % 0.0047) rn.- 1 . (24) 

The pion absorption may be seen as contributing roughly 10% to the real part. 

A couple of recent papers have reported an extraction of A n - d from the pionic deuterium atomic transitions, Refs. 
[33l l8ll 182^ , using the Deser-Trueman formula [2!| to find (Coulomb interaction included) : 

A%- d = (-0.0259 ± 0.0011) + i (0.0054 ± 0.0011), (25) 

from Ref. [Hill, and 

A c n - d = (-0.0261 ± 0.0005) + i (0.0063 ± 0.0007), (26) 

from Ref. ; 82j, both in units of the inverse pion mass as before. Both of these data are consistent with each other, 
and particularly the imaginary parts are also consistent with the model calculations in Refs. [78l l80j|. However, the 
real parts are about half the magnitude of the model result in Ref. 80] . If we assume that the experimental result be 
correct and that the result of the model calculation mentioned above be also correct regarding the sign and size of the 
absorption contribution to the real part, then the real part of the scattering length without the pion absorption effect 
should be about —0.020 m^ 1 . Then it appears that the data indicate the pion absorption effect to be very close to 
30%. And supposing that we translate this to our present K~d scatter ing leng th, the pure three-body result could not 
be acceptable. But there is a possible way out of this impasse: in Ref . |33|, |8l| , by combining the data from the pionic 
hydrogen and pionic deuterium, the value of the scattering length 7r~n, or more precisely the isoscalar combination: 

26§ = + (27) 

turned out to be consistent with zero. This near vanishing of &5 should be expected from current algebra calculations 
|83j |. and particularly in the soft pion (viz. zero pion mass) limit (however, we should be reminded in this respect 
that the extraction of this latter quantity is still rather model dependent, see Ref. [^3, US for details, as well as the 
consequence from the ttN partial wave analyses, see for example Ref. 84]). So if one accept this result, the lowest 
order pion-deuteron scattering length in a static calculation vanishes, and this is the basic origin of the smallness of 
TZe(A 7l - d ). Consequently, one might well come up with that large pion absorption effect. 

Now, this is far from true in the case of the K~d scattering length where (i) because the kaon cannot be regarded 
as soft, and (ii) because of the predominantly exo-energetic nature of the associated coupled channels, even the lowest 



order scattering length is neither vanishingly small nor purely real to start with, even without K~ absorption. Hence 
we may safely abide by the estimate of the kaon absorption effect as discussed earlier, and so we set the effect to be 
less than 10%. Of course a more quantitative study will have to be done. 

Our present calculation has not included any electromagnetic interactions, of which the Coulomb interaction plays 
the dominant role in the actual hadronic quantities measured. So we now come to discuss the effect of the Coulomb 
interaction as our last subject for this Section. Here just like what we have done above, we will borrow a good part of 
our argument below from the ir~p and n~d scattering length problems. In fact, as long as the aspect related to the 
electromagnetic interactions is concerned, replacing 7r~ by K~ should not alter it in an essential manner. So the first 
important point to be stressed is the following: all the experimental determination of the ir~p elastic (and ir~p — > ir°n 
charge exchange) scattering lengths extracted from the pionic hydrogen atom level shift and width have taken care of 
various electromagnetic corrections (including the finite electromagnetic size of the pion, vacuum polarisation effect, 
etc.) to the Deser-Trueman formula only [33L l8ll l82|. As found, for example, in Ref.[85j], this correction is up to 
about 2%. But the effect of the point Coulomb interaction, which is the very basis for the use of the Deser-Trueman 
formula, has not been taken out. Thus to be more precise, the extracted quantities should carry an index "c" to 
indicate that the Coulomb effect is still there. This we have done explicitly in what we have written above. In the case 
of the pion-deuteron scattering length, even this type of electromagnetic corrections to the Deser-Trueman formula 
has not been attempted. To a large extent the reason should be that the calculation is far more complicated than 
for the Tr~p case. But one may well suspect that a straightforward application of the Deser-Trueman formula to this 
already explicitly extended system would obviously introduce an error far larger than these sophisticated correction 
to that very formula. 

There have been model dependent but rather detailed calculations relating a c K - p and a w - p , the latter being due to 
the purely strong interaction [27, 86]. It was stated that the difference between the two quantities is just a fraction of 
a percent. As the same method cannot b e ap plied, a very simply estimate was carried out to relate the corresponding 
quantities for the tt~ d scattering length |27|: 

POO POO 

K-d/A~d= u 2 Deu (2r)<p 2 (0,r)dr/ u 2 Deu (2r)r 2 dr, (28) 
Jo Jo 

where U£> eu (r) is the S- wave radial deuteron wave function and 0o(O, r) is the S'-wave zero momentum Coulomb wave 
function for a unit charge. Clearly, this is just to semi-quantitatively introduce the distortion of the incoming charged 
pion due to Coulomb interaction. The result is about 4% increase in the magnitude for A v - d as compared with its 
Coulomb included counterpart. However, as we see in the experimentally extracted AZ,_, reported above the 
error bars are just about the size estimated here. So here again, the Coulomb correction seems to appear quite small. 
Likewise, the same line of reasoning might well apply to the case of the K~d scattering length. When translated into 
the model prediction, a possible allowance should be taken into consideration between the purely strong and Coulomb 
included scattering lengths, although by its very nature the estimate should be regarded qualitative. To this end it 
should be useful to refer to the work of Barrett and Deloff (34|. They introdued a set of rather simple K~d optical 
potentials and calculated the strong interaction shift and width of the IS atomic level for the kaonic deuterium. Also 
the optical potentials were used to calculate the purely strong (A K - d ), as well as the Coulomb included (A c K _ d ) K~d 
scattering lengths. The observation they made was that the Deser-Trueman formula might be inaccurate, and that a 
blind application of that formula and the identification of the extracted quantity as the K~d scattering length due 
only to strong interaction might introduce an error as large as 20%. A word of caution should be due regarding this 
work: the optical potentials constructed there were quite simple, and strong non-locality expected from the dominance 
of the A(1405) was absent. So a more realistic optical potential should be constructed in order to give more reliable 
statements on the issues. Otherwise, if one wants to simply obtain the Coulomb included K~d scattering length, it 
is possible to use the pure-strong K~d amplitude and apply Coulomb corrections as found, for example, in Ref. 87]. 



VI. SUMMARY AND CONCLUSIONS 



The scope of the present work being to develop a reliable formalism to calculate the K~d scattering length, our 
starting point was a thorough study of various appropriate two-body processes. For that purpose, we first focused on 
the highly inelastic K~p initiated reactions in the kinematics region p l £ b < 250 MeV/c. The elastic, as well as the 
relevant seven inelastic coupled-channels were investigated via an effective non-linear chiral meson-baryon Lagrangian. 
Within a broken SU(3)-symmetry scheme, the adjustable parameters of the formalism were determined by a fitting 
procedure on the threshold branching ratios and total cross-section data, leading to reduced x 2 's close to 1.2. 

To make clear the sensitivity of the observables to the phenomenological ingredients, three models were constructed. 
They were then exploited to predict other measured quantities, namely, 



• K p scattering length, for which our best value is 

a K - p = (-0.90 + i 0.87) fm, 
in agreement, within the experimental uncertainties, with the recent KEK data, 

a c K - p = (-0.78 ± 0.15 ± 0.03) -H(0.49± 0.25 ±0.12) fm. 

• The E-7T mass spectrum, measured at CERN some 20 years ago, was reproduced in an acceptable manner. This 
quantity was found quite sensitive to the rjY intermediate state within the used coupled-channel approach. 

In view of the K~ d system investigations, besides the K~p interactions, one of course needs another elementary 
amplitude: K~n, for which no data is available. We hence performed predictions for the elastic and the four 
inelastic coupled-channels. Moreover, predictions for various scattering lengths were made, i.e., a n (K~n — > K~n), 
a°(if n — > K n), and a ex (K~p — *■ K n). Finally, the implementation of the two-body sector was completed by 
studying other relevant channels: pion-nucleon, nucleon-nucleon (the deuteron), and nucleon-hyperon interactions. 

Then, we moved to the central topic of the present work and developed a relativistic version of the three-body 
Faddeev equations, to which we embodied the above elementary operators. As expected, this formalism allows us to 
go far beyond previous investigations, such as single- and multi-channel approaches and fixed center approximation. 

As in the case oiKN interactions, we performed our studies in both isospin and particle bases. 

Investigating the K~d system brings in phenomena with small contributions. A quantitative determination of 
their importance goes beyond the scope of the present work. However we tried to evaluate them qualitatively, i.e., 
possible contributions from four-body intermediate states, K~ absorption, and Coulomb correction, this latter being 
also present in the elementary two-body channels. 

Finally, our best value for the K~d scattering length is, 

A K - d = (-1.80 + i 1.55) fm. 

Given the quality of the phenomenological input and approximations introduced, our estimations lead us to attribute 
to the above values an uncertainty of about 10%. 

The awaited for data will soon make clear how realistic our predictions are. These experimental results will come 
from DA&NE on the KN and K~d scattering lengths, as well as from COSY, ELSA, JLab, and SPring-8 on the 
lowest mass A-resonances, including the T,n mass spectrum. From theoretical side, several topics deserve to be studied, 
such as Coulomb effects and K~d break-up channels. 
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APPENDIX A 

Particle and isospin basis 

We consider two particles with isospin Ii, I2, projections I± z , I2Z, and total isospin J, I z . The transformation from 
isospin to particle basis is written as: 

I1+I2 

\hhzhhz)= (hhzhhz\Ihz+hz)\Ihz+h Z ). (Al) 

I=\h-Ii\ 

Next, we specify the values of the isospins and their projections of the particles, following the phase convention 
given in Tabled for the isospin states. Calculating the appropriate Clebsch-Gordan coefficients, we easily obtain the 
linear relations between the states in the two bases. 



1. KN states 



For the K p interactions, the eight physical states defined in Eqs. JSJ) are expressed as linear combinations of the 
7 = 0,1 and 2 states, according to the following relations: 
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From these expressions, we can deduce the relations between the transition potentials in the two bases. Choosing 
a separable form as Eq. J2J), with isospin- independent form factors, we obtain the relations between the strength 
parameters in the two bases. Neglecting the contributions of the 1 = 2 states, we have for example: 



^K-p-K°n 2^ KN ~ KN ^KN-Kn\ ' \r°ri-7r+S- /g "ViV— jrE 2 ^KN—irS ' (^) 

and so on for the other parameters, with the symmetry property: Ay = Xji. 

In the case of the K~n and related states, only the 1=1 and 2 states contribute, and the five physical states Eq. 
© are expressed as follows: 



' \ K-n) = -|1 -1) FN , 

Itt-E ) = -^11-1)^ + ^12-1)^, 
< |7r°E-) = -^11-1)^ + ^12-1)^, 

I 7T-A) = |1 -l) wA , 

J»7S-) = |1 -1), E . 

The relations between the strength parameters in the two bases are obtained along the same lines as in the K~p 
case. Neglecting as before the 1 = 2 states contributions, we obtain for example: 
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KN-KN 
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(A3) 



etc. 



2. 7riV states 



The different ttN states in the particle basis are expressed as linear combinations of the J = i, | isospin states. 
Considering the possible charge states of the pion and nucleon, we have the following two sets of coupled states: 



and the following single state: 



k = I ! - |w 

Retaining only the contributions with total isospin 3/2, we deduce from the above expressions the following relations 
between the strength parameters in the two bases: 

2 1 

etc. 

3. EiV-ATV states 

The UN- AN states in the particle basis are expressed as linear combinations of the / = |, | states. We have two 
sets of three coupled states: 

o„\ _ . / 2 I 3 1 \ i 1 1 1 



and the two following single states: 



l S p } = V 3 I 2 2 ) SJV ~ 2)SJV, 

I E+n) = I f i) SA r - y§ I i |)ejv, 
|Ap) = 1 1 1 >AiV, 



An) = || - i) A 7v- 
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2 2 >CT ' 



From these relations, we can express the strength parameters in the particle basis in terms of those in the isospin 
basis. For example, we have: 

\/2r>l/2 _ \ 3 /2 1 \ _\ X 1 / 2 fAtn 

g L A £iV-£iV A EiV-£ivJ ' i A S+ri-Ap — A S-p-An — U jj A SJV-AJV' l > A0 /' 



and so on. 

APPENDIX B 

Separable model for coupled channels 

The Bethe-Salpeter equation |42j is the relativistic generalization of the Lippmann-Schwinger equation describing 
the scattering of two particles. In the case of n coupled two-body channels, we have a system of n coupled equations, 
which reads, in operator form : 

Tij(a) = Vij + Y, V lk G k a (a)T kl (a), (Bl) 

k 



with a the square of the total centre of mass energy. The indices {i, j, k} run over the n two-body channels. Go is the 
two-body propagator, Vij the transition potential between channel i and j, and T^- the i-matrix for that transition. 
After projecting in the 4-momentum space representation, we obtain a system of 4-dimensional integral equation. 



Using the Blankenbecler-Sugar (BbS) method (Refs. |4j, |4J, |45j ) , this equation can be reduced to the following set 
of coupled three-dimensional equations: 

Tij(pi,Pj;<r) = %(pi,Pj-;o-) +^2 d 3 p k Vik{pi,Pk] v)G%(p k ; a)T kj (p k ,pj; er), (B2) 

k J 

where pi is the momentum in channel i. The relativistic two-body propagator has the following relativistic expression: 

G »<^> - i^2^ra^r (B3) 



where the two particles in channel k are labeled as 1 and 2, and = y pf + ntf is the energy of particle i with mass 
mi. 

Finally, in the case of s-wave interactions, Eq. I|B2|1 reduces to the following set of coupled one-dimensional integral 
equations: 

poo 

T t3 (p t ,Pj]a) = I ,.,(/',. /',.:t; + ^2 Pk d PkV tk (Pt,Pk;<j)Go(p k ; a)Tkj{j>k,Pj]o-). (B4) 

k Ja 

Now, we assume that Vy's are separable, i.e. we write: 

Vij (pi , pj ; a) = gi (pi) Xij (a) g 3 ( Pj ) . (B5) 

The g's are the form factors, and is the strength for the transition i «-> j. Here we assume that in general A's 

are functions of a. Using this expression in Eq. (|B4|I . we obtain T as separable: 



Tij{pi,Pj-,a) = g l {p i )R l j{a)g j {p :j )i (B6) 
where Rij is an element of the following (n x n) matrix: 

R(a) = [X- 1 (a)-G(a)}-\ (B7) 
Here A(<r) is an n x n matrix of the strengths, and G(s) is a diagonal matrix with the following elements: 

/>oo 

Gij(a) = % / p 2 dpg 1 l (p) G*(p; a), (B8) 
Jo 

where Gq is the relativistic two-body propagator for channel i, calculated from Eq. (|B3(1 . 

Note that in the case where the strength matrix has no inverse, Eq. (|B7J| must be re-written as: 

R(a) = [1 - \(a)G(a)}~ 1 \(a), (B9) 



where 1 is the (n x n) unit matrix. This situation occurs for example when the KN interactions are considered in 
the particle basis, as the strength parameters obtained as linear combinations of the isospin basis values constitute a 
matrix with its determinant equal to zero. 



APPENDIX C 



Antisymmetrization 

In the isospin basis, the two nucleons are considered as identical particles. Therefore, one must construct Born 
terms and three-body amplitudes properly antisymmetrized with respect to the two nucleons. We present hereafter 
the principle of the method, and we refer the reader to Refs. 0, 0] for more details. 

Let us label the two nucleons as N± and N 2 . The d channel corresponds to K{N\N 2 ), where the tilde means that 
the deuteron wave function is properly antisymmetrized. We must also add the nucleons labels to the labels defined 
in Table IVTTl A* = T,(wNi), ft = n{YNi), a 1 = N % (ttY), with i = 1,2, and y l = N {KN % ), with i,j = 1,2 (i ^ j). 
(note that the lower indices have been removed, since they do not participate in this discussion). 

Now, we have two types of Born term, depending on whether one or two nucleons are involved. In the first category, 
we have Za , and Z a p. For these terms, we only need to introduce the nucleon index. For example, we define: 

Z A x a i = (£(7riVi) |Go| iVi(7rS)), Z A 2 a 2 = ( ZJ(7rJV 2 ) |Go| N 2 (wT,))., where G is the three-body propagator. Since 
nucleons 1 and 2 are identical, it is clear that these two Born terms are identical, and we will set: Z Aa = Z A i a i = 
Z A 2 a 2. In the same way, we define: Z A p — Z A tni — Z A 2^2, and: Z a p = Z a ipi = Z a 2^2. 

The Born terms involving two nucleons, namely Zd y and Z yy , must be antisymmetrized. As the d state is already 
antisymmetric, we only need to antisymmetrize the y state. So, using the notations defined above, the antisymmetric 
Zdy Born term will be defined as: 

Z dy = (K(Kn 2 ) \G q ^= [I N 2 (KN 1 ) ) - | Nx(KN 2 ) )] = ^[^V - z dyA- (CI) 
Exchanging N\ and N 2 in Zd y i, it is obvious that: Zd y i = —Zd y 2, thus: 

Zdy = V2Z dy i . (C2) 
We proceed along the same lines to define the antisymmetric Z yy Born term: 

Zyy = [{ N 2 (KN 1 ) \ - ( N 1 (KN 2 ) |] G ^= [| N^KN,) ) - | Nx(KN 2 ) )] . (C3) 
AsK can be exchanged only between the (KNi) and (KN 2 ) pairs, we have: 

Z yy = -\ [( N 2 {KNi) \G \ N^KN^ ) +(N 1 (KN 2 ) \G \ N 2 (KNi) )] = \\Z yly 2 + Z y2yl }. (C4) 
Due to the identity of the two nucleons, we have: Z y i y 2 = Z y 2 y i, and thus: 



Zyy — —Z y l y 2. (C5) 

For the practical calculation, only the coefficients appearing in the expressions of the antisymmetric Born terms 
are important, and we can ignore the nucleons labels. 

Concerning the two-body propagators, the nucleons labels can be ignored, since only one nucleon is eventually 
involved in the propagating pair, (the deuteron propagator corresponds to a properly symmetrized 3 Si or 3 Si- 3 Di 
state). 

Finally, the three-body equations can be rewritten in the same form as Eq. l|13Jl . where now the Born terms are 
the antisymmetric terms as defined above. 
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